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Abstract 

A simple fast algorithm of the conodes calculation in binary sys- 
tems is proposed. The method is based on exact solution of the prob- 
lem on common tangent to pair of approximating parabolas. Sequence 
of approximating parabolas pairs having second order tangency points 
on separated concave parts of free energy isotherm is generated. The 
sequences of the tangency points pairs conodes with approximating 
parabolas converge to equilibrium compositions of phases in binary 
systems. 



1 Introduction 

Free energy of a binary system at fixed volume V or pressure P and 
temperature T depends on the system composition x only. A typical form of 
function graph in two-phase region presented on Figure [T] Equilibrium 
compositions xi, x<i of the coexisting phases in binary systems by known 
free energy are determined by tangency points of the cono d<£] with 

function graph. 

Thus, the problem of x,\ and Xi evaluation reduces to pure mathematical 
task of the conode for given function finding. There are several methods 
of this task solutions [2J. Most of them require a number iterations to have 
accuracy of the order of 1CT 3 . 

1 conode is the common tangent to separated concave parts of function [TJ p. 71] 
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Figure 1: Schematic free energy — composition x diagram for a binary 
system in two-phase region. 
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The aim of this paper is development of an algorithm for evaluation of the 
equilibrium compositions of phases in binary systems with same accuracy in 
2-3 iterations. 



2 The algorithm 

The algorithm of the problem solution consists the following steps 

1. Choice in the zero-order approximation some arbitrary points 

in first and second concavity intervals, respectively. These points are 
the second order point of tangency with tangent parabolas. 

2. Derivations of the equations of tangent parabolas to curve $(x) in 
points x^\ respectively. 

3. Exact analytical finding of the conode to approximating parabolas and 
determination the point of contact conode with approximating parabo- 
las. 

4. Choice in the first- order approximation x-coordinates I, of the 
point of contact conode with approximating parabolas and transfer to 
point. 

Performance of iterations stops after achievement of the set accuracy. 



3 Analytical calculations 

In zero approach we shall choose points on the different parties of a convex 
site of a curve <&(#). The equations of parabolas approximating free energy 
in vicinities of points xf\ (i — 1, 2) have the following form 

fW(z) = a? + b?(x-x?) + C ^(x-x?y ; (1) 

(0) 9 

/««» M = a?> + b?> (x - 4°») +\(*- 4°») , (2) 

4°> = *«h "f = *'(4 0) ); «£" = •"(*?) (3) 

are values of function $ (x) and its derivatives of the first and second orders 
in points x^ and xf\ The equations of tangents to each of parabolas in 



where 
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corresponding points x± and x? have the following form: 
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The condition of coincidence of these tangents corresponds to system of 
two equations concerning two abscissas points of a contact x? and x 2 to 
the general tangent to both approximating parabolas Q and ([2]): 
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The solution of system (J6J) in terms of a deviation from initial approxi- 



mation Ax 
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x- ) leads to following result 
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and Ax? satisfies to the quadratic equation 

A«(Axi 1) ) 2 + J B«Ax« + C'« = 0, 
with following coefficients 
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It is necessary to choose from two solutions of the equation (|8j) the solution 
that has a finite limit at {c\ — c 2 ) — > 0: 
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The points x[ and x 2 finding indicates the first iteration ending. In the 
second iteration these points x± and x^p play the same role as the points x± 
and in the first iteration and so on. 

There are a number ways or the algorithm accuracy setting. The simplest 
way of the iteration termination is stopping of the process by the following 
condition realization 

< e, (13) 

where e is a some preassigned value. In particular, it can be precision of 
experimental data or any other quantity. 

Let us consider some simple examples of the offered algorithm realization. 

4 Application of a method to some forms of 
the thermodynamic potential 

As examples, we consider the following model expressions for free energy: 

1. Free energy in a polynomial of the forth power form as a simplest 
example; 

2. Van der Waals free energy; 

3. Redlich-Kister form of free energy. 

Realization of the algorithm fulfilled in Mathematica 5.0 system, results out- 
put for graphical representation carried out in Excel system. 
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4.1 A polynomial form of thermodynamic potential 

The simplest model potential having characteristic for a free energy in two- 
phase region behavior is the polynomial of the fourth power of x. This aim 
was achieved by the polynomial coefficients selection. Let us consider one of 
such polynomials 

$(x) = -x 4 - — x 3 + x 2 - -x. (14) 

Graph of this function is presented on figure ([2]). 

The calculations show, that the relative error in determination of the 

I (3)_- I 

tangency points after three steps does not exceed value — - < 5 ■ 1CT 4 , 
and after fourth step the error is less as 1CT 7 independent of the initial 
points xf \ x 2 choice. 
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Figure 2: Search of the conode to model free energy in the form of the fourth 
power polynomial = \x A — |^x 3 + x 2 — |x. Notations: the black line is 
the function graph; colored dotted lines are graphs of the approximating 
parabolas in different approximations; straight lines are the tangents; the red 

(k) (k) 

points are the points x\ , x 2 , the pricked red points are the points of contact 
between common tangents and approximating parabolas; the yellow points 
are the points of contact between conode and function 
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4.2 The Van der Waals free energy 

The Van der Waals free energy has the following form: 

G(P,T,N) = —RTN ln(V — Nb) - + PV + c(T, N), (15) 

where a, b are the Van der Waals constants, c(T, N) is an independent of 
the volume function. Let us introduce a new variable x = Nb/V: 

_ rl {Nb(l-x)\ aNx bPN 

G(x) = —RTN In i '- — + + c{T,N) 

V x J b x (lg) 

f 1 — x\ a 3 v ' 

= ${x) = — ai In — a<ix H h 04, 

\ X / X 

where a« are the constants related with the Van der Waals parameters and 
variables N, T, P. Results of the calculations for case a± = 2019.8, 02 = 
8531, 0,3 = 85.34, = 21323 (these parameters correspond to carbon dioxide 
C0 2 at temperature T ~ 0.8 T c and pressure P ~ 20 bar) presented on 
figure ([3}. 

4.3 The Redlich-Kister potential 

The Redlich-Kister potential 

n 

$ = ^2x l $ i + RT^Xilnxi + ^XiXjLij^ixi- Xj) k (17) 

i i j>i k=0 

is one of the most often used for modeling of the multicomponent systems 
thermodynamic properties. 

The calculations results for the Redlich-Kister potential 

$(x) = 2.5 [xlnx + (1 — x) m(l - x)] 

(18) 

+2.2x(l - x) (1 + 2x) + 0.2x + 0.8(1 - x) 
are presented on figure Q. 

5 Conclusion 

The presented method for determination of coexisting phases compositions 
qualitatively differs from usually used methods based on a free energy nu- 
merical minimization [21 [3J HI |5j EJ El [9] . 
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Figure 3: The calculations results for Van der Waals free energy $(a;) = 
-2019.8 In ((1 - x)/x) - 8531a; + 85.34/x + 21323. The notations are the 
same as on figure |2j 
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Figure 4: Search results of the conode for Redlich-Kister potential $ (x) = 
2.5 [x \nx + (1 - x) ln(l - x)} + 2.2x(l - x) (1 + 2x) + 0.2x + 0.8(1 - x). The 
notations are the same as on figure [2} 

Using this method to several model thermodynamic potentials has shown 
the efficiency of the algorithm. As a rule, it is enough to three iteration to 
obtain accuracy of the coexisting phases compositions up to 10 -4 . Each of 
the iterations includes nine computing operations (evaluation function and its 
first and second derivatives in the initial points n , and then Ax' ra+ > , 

x^ n+1 ^ X2™ +1 ^). The rigorous analytical estimate of degree of convergence will 
be derived later. 
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